Goodness-of-fit tests for neural population models: the multivariate time-rescaling theorem

Poster Presentation from Nineteenth Annual Computational Neuroscience Meeting: CNS*2010 San Antonio, TX, USA. 24-30 July 2010 Statistical models of neural activity are at the core of the field of modern computational neuroscience. The activity of single neurons has been modeled to successfully explain dependencies of neural dynamics to its own spiking history, to external stimuli or other covariates [1]. Recently, there has been a growing interest in modeling spiking activity of a population of simultaneously recorded neurons to study the effects of correlations and functional connectivity on neural information processing (existing models include generalized linear models [2,3] or maximum-entropy approaches [4]). For point-process-based models of single neurons, the time-rescaling theorem has proven to be a useful toolbox to assess goodness-of-fit. In its univariate form, the time-rescaling theorem states that if the conditional intensity function of a point process is known, then its inter-spike intervals can be transformed or “rescaled” so that they are independent and exponentially distributed [5]. However, the theorem in its original form lacks sensitivity to detect even strong dependencies between neurons. Here, we present how the theorem can be extended to be applied to neural population models and we provide a step-by-step procedure to perform the statistical tests. We then apply both the univariate and multivariate tests to simplified toy models, but also to more complicated many-neuron models and to neuronal populations recorded in V1 of awake monkey during natural scenes stimulation. We demonstrate that important features of the population activity can only be detected using the multivariate extension of the test. ..

Incremental Mutual Information: A New Method for Characterizing the Strength and Dynamics of Connections in Neuronal Circuits

Understanding the computations performed by neuronal circuits requires characterizing the strength and dynamics of the connections between individual neurons. This characterization is typically achieved by measuring the correlation in the activity of two neurons. We have developed a new measure for studying connectivity in neuronal circuits based on information theory, the incremental mutual information (IMI). By conditioning out the temporal dependencies in the responses of individual neurons before measuring the dependency between them, IMI improves on standard correlation-based measures in several important ways: 1) it has the potential to disambiguate statistical dependencies that reflect the connection between neurons from those caused by other sources (e. g. shared inputs or intrinsic cellular or network mechanisms) provided that the dependencies have appropriate timescales, 2) for the study of early sensory systems, it does not require responses to repeated trials of identical stimulation, and 3) it does not assume that the connection between neurons is linear. We describe the theory and implementation of IMI in detail and demonstrate its utility on experimental recordings from the primate visual system

Single-neuron dynamics in human focal epilepsy

Epileptic seizures are traditionally characterized as the ultimate expression of monolithic, hypersynchronous neuronal activity arising from unbalanced runaway excitation. Here we report the first examination of spike train patterns in large ensembles of single neurons during seizures in persons with epilepsy. Contrary to the traditional view, neuronal spiking activity during seizure initiation and spread was highly heterogeneous, not hypersynchronous, suggesting complex interactions among different neuronal groups even at the spatial scale of small cortical patches. In contrast to earlier stages, seizure termination is a nearly homogenous phenomenon followed by an almost complete cessation of spiking across recorded neuronal ensembles. Notably, even neurons outside the region of seizure onset showed significant changes in activity minutes before the seizure. These findings suggest a revision of current thinking about seizure mechanisms and point to the possibility of seizure prevention based on spiking activity in neocortical neurons

Exponential martingales and changes of measure for counting processes

We give sufficient criteria for the Dol\'eans-Dade exponential of a stochastic integral with respect to a counting process local martingale to be a true martingale. The criteria are adapted particularly to the case of counting processes and are sufficiently weak to be useful and verifiable, as we illustrate by several examples. In particular, the criteria allow for the construction of for example nonexplosive Hawkes processes as well as counting processes with stochastic intensities depending on diffusion processes

Intrinsic gain modulation and adaptive neural coding

In many cases, the computation of a neural system can be reduced to a receptive field, or a set of linear filters, and a thresholding function, or gain curve, which determines the firing probability; this is known as a linear/nonlinear model. In some forms of sensory adaptation, these linear filters and gain curve adjust very rapidly to changes in the variance of a randomly varying driving input. An apparently similar but previously unrelated issue is the observation of gain control by background noise in cortical neurons: the slope of the firing rate vs current (f-I) curve changes with the variance of background random input. Here, we show a direct correspondence between these two observations by relating variance-dependent changes in the gain of f-I curves to characteristics of the changing empirical linear/nonlinear model obtained by sampling. In the case that the underlying system is fixed, we derive relationships relating the change of the gain with respect to both mean and variance with the receptive fields derived from reverse correlation on a white noise stimulus. Using two conductance-based model neurons that display distinct gain modulation properties through a simple change in parameters, we show that coding properties of both these models quantitatively satisfy the predicted relationships. Our results describe how both variance-dependent gain modulation and adaptive neural computation result from intrinsic nonlinearity.Comment: 24 pages, 4 figures, 1 supporting informatio

Human seizures couple across spatial scales through travelling wave dynamics

Epilepsy—the propensity toward recurrent, unprovoked seizures—is a devastating disease affecting 65 million people worldwide. Understanding and treating this disease remains a challenge, as seizures manifest through mechanisms and features that span spatial and temporal scales. Here we address this challenge through the analysis and modelling of human brain voltage activity recorded simultaneously across microscopic and macroscopic spatial scales. We show that during seizure large-scale neural populations spanning centimetres of cortex coordinate with small neural groups spanning cortical columns, and provide evidence that rapidly propagating waves of activity underlie this increased inter-scale coupling. We develop a corresponding computational model to propose specific mechanisms—namely, the effects of an increased extracellular potassium concentration diffusing in space—that support the observed spatiotemporal dynamics. Understanding the multi-scale, spatiotemporal dynamics of human seizures—and connecting these dynamics to specific biological mechanisms—promises new insights to treat this devastating disease

Differential Cytokine Gene Expression According to Outcome in a Hamster Model of Leptospirosis

Leptospirosis is a widespread bacterial infection that is transmitted by soil or water contaminated by the urine of infected animals, or directly from these animals. It has highly diverse clinical presentations, making its differential diagnosis difficult. Though most cases are minor and self-resolving, there are also severe forms that include a sepsis pattern and multiple organ failure, and have possible fatal outcomes. Predictors of disease evolution and outcome are scarce, yet they would be very valuable to clinicians as well as to better decipher disease pathogenesis. In this study, we used a hamster model of leptospirosis to evaluate if immune genes were differentially expressed between individuals and if their expression levels could help forecast the outcome of the disease. We found that hamsters that later died from leptospirosis had significantly higher expression levels of both pro- and anti-inflammatory mediators compared to survivors. These results suggest that expression levels of these immune effectors might be helpful predictors of outcome in leptospirosis and that septic shock contributes to fatal leptospirosis